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ABSTRACT 



The spectrum of the stationary synchronous interval process in 



the stochastic point process obtained by superposing p Erlang renewal 
processes is derived by using relationships based on the Palm-Khinchine 
formulae and the fundamental identity linking the counting process of a 
point process to the interval process. The spectra coincide with those 
of mixed moving average — autoregressive processes. Explicit results are 
derived for a few simple cases for small p and a computational formula 
for the more complicated cases. Some general results on the shape of 
the spectrum of intervals are also given. 



Prepared by: 



I. Introduction. 



Much work has been done recently in the area of the superposition 
of independent stationary point processes, extending earlier results 
summarized in Cox and Lewis (1966, Ch. 8) and Cinlar (1972). Various 
approaches have been utilized in deriving the probabilistic properties 
of the superposition process, including the development of expressions 
for the description of counts of events occurring in some interval and 
for the description of the properties of intervals between events 
starting at an "arbitrary” event. Following Lawrence (1971) we call 
these the synchronous intervals of the process. 

The present paper considers the problem of describing the second- 
order joint moment structure of the synchronous intervals for the special 
case of the superposition of identical Erlang renewal processes. 

Notationally , let p be the number of component processes to 
be superposed and let each component process be an Erlang renewal pro- 
cess, i.e. have independently distributed intervals between events, each 
interval having the same Erlang probability density function 



The superposition of p of these independent processes gives the super- 
posed process to be studied. We use the description "Erlang process" 
rather than "Gamma process" to indicate that the shape parameter k 
in (1.1) is an integer. 




k-1 -kx/y 
X e 



(x ^ 0). 



( 1 . 1 ) 



r(k) 



Downton (1972) considered the problem of deriving the distribu- 
tions of synchronous intervals in this superposition process and utilized 
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both an imbedded Markov property and an integral equation approach to 
develop limiting forms for the distributions, i,e. for large times of 
observation and large numbers of superpositions, Downton was concerned 
with applications to reliability studies; note that the superposed 
process is the overall failure pattern for a system consisting of p 
renewable components operating in series, i.e. the system fails whenever 
a component fails and is renewed. 

Lawrance (1973), amongst other things, extended Downtons results 
and showed that the synchronous intervals between events in the super- 
posed process are correlated. He calculated the first two serial 
correlation coefficients for specific cases of superposed Erlang renewal 
processes, and obtained joint and univariate distributions of synchronous 
intervals . 

Barnett (1970) used simulation to investigate some of the 
correlational properties of superposed processes in the specific area 
of "studying the supply pattern of goods delivered to a supply depot 
by a fleet of vehicles." His processes may be considered to be Erlang 
processes with large values of k. 

Here we develop somewhat more general expressions than those 
given previously for the serial correlations and the spectrum of syn- 
chronous intervals in superposed Erlang renewal processes. Specific 
forms for the spectra are given for a few cases. A general expression 
for the generating function of the distribution of counts is used to compute 
the spectrum of intervals and serial correlations in other cases. Some general 
conclusions about the shape of the spectrum of the intervals in the super- 
posed process are given; these utilize results of Enns (1970) to the 
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effect that if intervals in the component processes have increasing 
failure rates (decreasing failure rates) , then so do the marginal syn- 
chronous intervals in the superposed process, Enns (1970) also discussed 
serial correlations and tried to relate it to the properties of the 
failure rates of the intervals in the component processes. 
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2 . Preliminaries . 

Cox and Lewis (1966, Ch. 4) used the fundamental identity of 
point processes, 



(N(t) < k) (S > t) (k = 1,2,...) (2.1) 

and Palm-Khinchine formulae to relate the generating function of the 
asynchronous counts, cf)(5 5 t) = to the spectrum of synchronous 

intervals {X^} in a stationary point process. The development leans 
heavily on Laplace transforms; for (1.1) we have 



g(s) = 



g(x)e dx 



= r_k/jj_'| 

'•s+k/y-' 



(Re(s) 2; 0) 



0 - 



= (^ S + l)-\ 



( 2 . 2 ) 



which is a rational function in s. 

Note that the parametrization for (2.2) gives an expected value 
of y; the intervals X can be construed to be the sum of k independ- 
ent exponentially distributed variables, each with mean y/k = g. 

Denoting the Laplace transform of the generating function for 
counts in an interval (0,t] following an arbitrary event (synchronous 
case) in the point process by 4>^(5;s), we have for a renewal process 
(Cox, 1962, p. 37) 






(C;s) 



1 - g*(s) 
s[l-{g*(s)}] 



(2.3) 



This expression is simple looking but in general difficult to invert. 
Using c()^(5,s) the Palm-Khinchine formulae (Cox and Lewis, 1966, Ch. 4) 
give the Laplace transform (£;s) of the generating function for the 
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number of events in (x,x-bt), where x is an arbitrary time in a 
stationary process (i.e. the asynchronous case). The relationship is 
(Cox and Lewis, 1966, Ch. 4, p. 68) 

(|)*(C;s) = ^ (|)*(5;s) . (2.4) 

If p renewal processes are superposed, we need to develop a general 

expression for the transform of the generating function of the asyn- 

* (p) 

chronous counts in the superposed process, denoted by <|) (£;s). 

This is because the spectrum can be obtained as (Cox and Lewis, 1966, 
p. 77) 

oo 

f . (co) = ^ {1 + 2 ^ p cos(Iod)} 

+ k=l ^ 

1 /(e^^;0+) + 

ir 2(()*(0,0+) - y ’ ^ 



where the are the serial correlation coefficients for synchronous 

intervals k lags apart, i.e. corr{X^,X^^^} , and the denominator is 
var(X^) /E(X^) = var(X^)/y. 

In principle it is easy, in the case of superposition of p 

*(d) *(p) 

identical processes, to obtain (J> ^ (5;s) or ^ (5;t)> since in 

the time domain 



= E{c ^ 



N, (t)+. . .+N (t) 

P } 



P N (t) 

= n E{? ^ } 

i=l 






( 2 . 6 ) 
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However, if <() (5;s) is known, as for the renewal process, to obtain 

^ (p) 

<t> (5;s) it is necessary to perform the operations (using 2.6) 

(p (S;s) -> (j)(g;t) (j) (5;t) ■> (p (S;s), 

where L ^ denotes inverse Laplace transformation and L a Laplace 
transform. Obviously the procedure would be cumbersome for large p, 
even if (p (5;s) could be inverted to give (j)(5;t). 

Section 4 of this paper is concerned with the development of 
<|)^^\5;t) or its transform in the special case of the superposition 
of independent, identical Erlang renewal processes, while Section 5 
develops the spectral density, for the superposed process. 

Section 6 gives specific functional forms for f^(w) and the 
p^’s in the cases of k = 2 and p = 2,3,4; the results are surpris- 
ingly simple. In Section 7 computational results for f^(^) are given, 
using the general results of Section 5, for k = 2,3 and p = 2, 3,..., 9. 
These illustrate the general form of the spectrum of the superposed pro- 
cess. Some results for superpositions of non-identical Erlang processes 
are also given, corresponding to cases discussed by Lawrence (1973) . 

In Section 3 we review some results of Enns (1970) and show how 
they can be used to derive some properties of the spectrum of superposed 
renewal processes, not necessarily Erlang. All the results are used in 
Section 8 in the discussion of the general shape of the spectrum of 
intervals of superposed Erlang renewal processes. 
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3 . Enn*s results and the initial point on the spectrum . 

It is clear that the marginal time between events in a stationary 
superposed process will have mean 

E(X^P^ = E(X)/p, 



where E(X) is the common mean time between events in the p component 
processes. Let the spectral density (2.5) of the synchronous intervals 
in the superposition process be denoted by (oj) and the coefficient 

of variation squared by = var (X^^^ ) /{E (X^^^ }^ . Then we have 

the following relationship between the asymptotic slope of the variance 
time curve (t) = var{N^^^(t)} and the initial point on the 

spectrum (Cox and Lewis, 1966, p. 78): 



lim 



(t) 

dt 



c^(x^^b 

E(X^^b 



(0+) . 



(3.1) 



Now = pV(t), where V(t) is the variance time curve for the 

component processes, and if the component processes are renewal pro- 
cesses it is well known (see Cox and Lewis, 1966, p. 81) that 



lim V (t) 



lim 

t-x» 



dV(t) ^ C^(X) 
dt E(X) ' 



Putting these results together we get, for superposed renewal processes, 




V («)E(X^Pb 

IT C (X^P^) 



p C^(X) ^ E(X) ^ C2(X) 

p c2(x^P^) IT c^Cx^^b 



There are two consequences of this result. 
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Enns (1970) showed that if X has an increasing hazard rate 
distribution (IHR) , then so does , and if X has a decreasing 

failure rate distribution (DFR) , so does X^^\ Thus, using results 
of Barlow and Proschan (1964) we have that 



Consequently 



c2(X^Ph if X is IHR, 



if X is DFR. 



ir C^Cx^^h 
1 

^2/v(p)' 



f_^Pho+) ^ 93S}Q . if X is IHR, 

+ 7T 



^ f^^hcH-) :S if X is DHR. 

+ 7T 



ir C^(X'‘^0 

Note that for the Erlang density function (1.1), X is IHR if k > 1 
and DFR if k < 1, and C^(X) = When k = 1 we have an exponential 
density and a Poisson process for the component processes and the super- 
position process; the spectrum then has the constant value 1 /tt and 
f^(0+) = fj^hcH-) = c2(x)/ir = l/ir. 

The second consequence of (3.2) is that in the limit as p 
when the superposed process goes to a Poisson process and C (X^^^) 1, 

the spectrum of intervals is tied down at o) = 0: 



limit f_^^^ (0+) 

p->oo 



C^(X) 

TT 



For the Erlang process this value is l/(Trk). Thus the spectrum, which 
is converging to the flat spectrum of the limiting Poisson process, has 
a ”notch" in it at o) near zero. This notch corresponds to the fact 
that if we could observe a superposed process well past the mean of 
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intervals in the component processes, the variance time curve would not 
be equal to t/E(X) , i.e. would depart from linearity. In other words 
low frequency effects and long time effects correspond. 

It is not known if f_^^^ (0+) goes monotonically from I/tt at 
p = 1 to l/(TTk) as p Computational results are given in 

Section 7. It would also be interesting to know the derivative of 
f_^^\u)) at 0 ) = 0+. 
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4 . General results for the generating function for an Erlang renewal 
process . 

For the Erlang renewal process, using (2.2) and (2.3) and letting 
X = k/p = 1/3, we have 

<('?(C;s) = r 

su - C(^» 



(X+s)^ - X^ 

Ic Ic ^ 

s{(X+s)‘" - X*"} 



(4.1) 



Hence, using (2.4), 



<{.*(5;s) = 



1 X(S-l) 
s ks 



(X+s)^ - X^ 
s{(X+s)'^ - ^x'^} 



(r){a+s)^ - + {^-^|^H(x+s)^ - 



a") 



(x+s)^ - 



(4.2) 



P(s) 
Q(s) ’ 



(4.3) 



where P(s) and Q(s) are polynomials in s. 

By expanding and collecting terms in P(s) one finds that 

P(s) = i Q(s) + (Q(s) + X^(C-l)} (4.4) 

, _i_ 2 _i_ _i_ k-1 

= a^ + a^s + + . . . + aj^_2S + aj^_^s , 



so that P(s) is a proper polynomial in s of order k-1, with 



a^ = X^"^{k + (C-l)(^)}, 
a^ = X^"^{^Hk+ (C-l)(^)}, 
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= + ( 5 - 1 ) 



\-2 "'' ^k ^ ’ 



Vl = 1* 



The first derivative of Q(s) with respect to s is 



Q'(s) = k(X+s)^ 



(4.5) 



Now let be the k distinct, complex-valued roots of 

5. i-e., 

Tj — 5 (j— 1,..* ,k) • 

The roots s^ of Q(s), i.e. clearly 






(4.6) 



and these too, like , are distinct and complex valued. 

We are now in a position to use a standard inversion formula 
(Gardner and Barnes, 1942, p. 338), for the ratio of two polynomials 
when the denominator has distinct roots: 

k P (s . ) s . t 



= y :iu_e j . 



We note that from (4.4) and (4.6) 



r-i 

P(s.) = (^^)2 (j=l,...,k), 

^ j 



and from (4.5) 



k-1 kX^ 

Q'(s.) = k(X+s.) = — (j = l,...,k). 

^ ^ j 



(4.7) 
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Consequently, the probability generating function for the asynchronous 
number of counts, N(t), in an Erlang renewal process is 



4>(C;t) = L ^{(j)*(C;s)} 



k P (s . ) s . t 

= y J e ^ 

A Q*(-) 









y — ^ 

3=1 



T, X(Xj-l)t 



7 e 



(4.8) 



We note too that the expansion for the variance-time curve of an 
Erlang renewal process obtained by Serfling (1970) is obtainable from 
(4.8). 
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5 . General Results for the Spectrum of Counts for Pooled Independent 
Erlang Renewal Processes . 

We consider here the case in which each of the p component 
Erlang renewal process have the same scale parameter X and shape 
parameter k. The probability generating function for counts in the 
pooled process starting at an arbitrary time is, from (2.6), 



Let J 
Section 4, 



be an index set. Using the result (4.8) from 








T. X(T.-l)t 



p 







T . 






e 



x{ I (t -1) }t 
j€J ^ 



Taking the Laplace transform of 






we have 



(5.1) 






(C-1) 



2 P 



k 

I 



k 

I 



T . 

3 



16J 

I S - * I (t -1) 

P ' j€J ^ 



(5.2) 



We note the following in (2.5): 



(i) 



/•7T 

f (o3)dto A 1, 

Jo 



(ii) the denominator in the second expression of (2.5) for f_^(m) 
is constant with respect to w . 
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Consequently, if we can evaluate (ff ^^\e^^,0+) + <p ^,0+) for 

0 0 ) ^ 7T, and integrate it, then (i) will give us the value of 

i7{2<fi *'^‘^(0,0+) - E(X)}. Thus we will have (o)) for 0 ^ w ^ ir, 

and also var(X^^^). Now when s = CH-, 






l)(C-l)2 P 

A n^( 



k 

I 

Jl=^ 



n 



MjkiiZr' 



(5.3) 



We are interested in the cases 5 = e^ and C = e We note that 



(€-1)^ ^ 

5 



5 - 2 + 5 



(5.4) 



Thus , 



(s-l)^ 



5=e 



^ (g-1)^ 

io) g 



g=e 



lU) o , -ItO 0/1 s 

. = e - 2 + e = -2(l-cos cu) 

- 10 ) 



notation 



1 

! _ 



1 th 
a k 


root of 


, th 
a k 


- io) 

root of e , 


, th 
a k 


- -i(jL) 

root of e 



In particular we assign indices to these roots in the following manner 

j = 1,. . . ,k, (5.5) 

J— l>...>k. (5.6) 



Xj = e.pu(^)). 



1,. =exp(-i(2l^)). 



Thus , ' 



T . 

J 



-1 
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Taking real and imaginary parts, 

Re(x.) = Re(T.') = cos , 

1 J ^ 

I(T.) = -I(Tj') = sin(^^^) . 

Because x.' = t. we see that 

1 J 

lU- = i ^ . 

(T.-l) T. - 2 + T.-l 



(5.8) 

(5.9) 



Thus, for any index set J, 



n 



= n 



... ( t .- 1)2 (x,'-l)2 

J €J J J 6 J J 



j6J e. - 2 + X. 
J J 



-L 



We now consider the terms 
in (5.3) when ^ is replaced by 



n . 

j€J -2{ 1-cos 

l/J^(l-x^) and l/^(l-x^') 

im , -im . ^ 

e and e respectively. 



(i) 



1 ^ 1 

I (1-x.) I {l-Re(x )}-i y I(t ) 

j6J ^ j6J ^ jSJ 

I U-Re(x )}+ i I I(x ) 

= j€J ^ j6J t . 

I I I U-Re(x )} - i I I(t.) I p 

j6J ^ j€J ^ 



(5.11) 



(ii) 



I (1-x!) I U-Re(x/)}-i I I(x ’) 

j6J ^ j6J ^ j6J ^ 

I U-Re(x .)} + i y I(t . ') 

^ i€J ^ j6J ^ 

II I U-Re(x ')} - i I I(t ')|P 

j€J ^ j6J ^ 
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I U-Re(T )}-i I I(t ) 
J€J ^ i€J ^ 

I I {1-Re(x.)} -i I I(t ) I I' 
j6J j€J ^ 



(5.12) 



Thus, for any index set J, 



1 (l-T.) I (1-t/) 



2a. 



j6J 



j€J 



(5.13) 



where 



and 



a, = ^ {1-Re(x.) } 

j6J J 

- I U-cos(2lf^)), 
j€J 

b. = I Kt.) 

^ j€J ^ 

= I . 

j6J ^ 



Let 



c = n {i-cos 
J6J 






(5.14) 



Gathering results (5.3), (5.4), (5.10), (5.13) and (5.14) we have 



fjP^oj) ~ (t)*^P\e^“,0+) + (t>*^P^(e ^“,0+) 



— ^ (1-cos 0 ))^ I 

Xk^P j^=l 



2a , 



j ^1 



(5.15) 



Because a^, and depend only on the indices in J and 

not on their order, the number of distinct terms to be computed for each 
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value of u) is equal to > the figurate number for k types 

of objects taken in sets of size p. The computation time for a , 

J 

hj and is proportional to p, so total computation time should 

be approximately proportional to p times the figurate number. 

A FORTRAN program was written to compute intervals 

of 0.05tt. Various values of k and p having figurate numbers no 
greater than 56 were used. Computing times were less than 2 seconds 
per run on the IBM 360/67 computer. Some results of the computations 
will be shown in Section 7. First we obtain some specific results 
for small p and k. 
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6 . Specific results for k = 2 and p = 2,3,4 , 

The generating function for the superposition of p 

identical Erlang renewal processes is analytically tractable for k = 2 
and small p. 

From (4.8), with a = and b = we 

obtain 

(KC;t) = (6.1) 

4C ' 

By taking {((>(S;t)}^ we obtain the generating functions for the 
superpositions of p such Erlang renewal processes. 

(i) For the case p = 2, we have 

{<l>(C;t)}2 = 4>^^^(c,t) = 

+ (l+C ^'^^)2 - 2 ( 1 - 0 ^ ( 6 . 2 ) 



whose Laplace transform is 



*( 2 ) 

<P 



(C;s) 



1 Hi-c^^^)* ^, (i+g^^^)** _ 2(1-^)^ 

I 65 I s - 2a s - 2b s-a-b 



(6.3) 



From (6.3) we can obtain the spectrum of intervals using (2.5). The 
result is 



fi^\oj) = - (1 - (0 ^ OJ ^ ir). (6.4) 

From (2.5) we have immediately that = -1/10, = 0 for k > 1, 

and the spectrum is that of a first order moving average process. This 
generalizes results of Lawrance (1973), who gave p^ = -1/10, but no 



higher correlations. 
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(ii) For p — 3 the generating function of the counting function of 
the superposition process is 

{_(i_5l/2)6 ^ 3 (i_^1/2)2(i_ 5)2 e(2a+b)t 

(4C ' ) 

- 3(1+C^^^)2(i- 5)2 e(a+2b)t (6.5) 



and its transform is 



s 3a s 2a b 



_ 3(l+C^^^)2(l-g)2 (i+g^/^)6 

s - a - 2b s - 3a 



.1/2, 



( 6 . 6 ) 



Again using (2.5) we obtain for the spectrum of intervals 

..(3)/ , 9 )121 - 48 cos 0 ) - 9 cos^ o)( , 

O (<»J) = n / (0 s: 0) s: tt) 

+ 25ir ) 41-9 cos m t 



(6.7) 



This is the spectrum of a mixed second-order moving average and 
first-order autoregressive process for which the serial correlations 
are expressed by three constants, since for k ^ 2 P,,-, ~ P3, • 

K. I _L K. 

Ic 

Equivalently, for k ^ 2, = C3 • 

In the present case = -0.1044447, ~ -.0316033, 

p^ = -0.0035118, or C = -2.55986 and 3 = 1/9. 



(iii) Skipping the details, we get for p = 4 the result for the 
spectrum 






8 J 


(s71 


- 502 cos w 


- 57 cos^ m - 24 cos^ o) 


379TTj 


1 


17 


- 8 cos 0) 



(0 ^ 0) ^ tt) (6.8) 
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This is the spectrum of mixed third-order moving average and first- 
order autoregressive process, so that the serial correlations may be 
expressed by four constants, since for k > 3, ~ Equivalently, 

for k ^ 3, 

For the case p = 4, k = 2 we get 



= -0,0979540 ; 

= -0.0442821 ; 

p^ = -0.0150263 ; C = -.961683 ; 
p^ = -0.0037597 ; 3=1/4. 

In Figure 1 we show the spectra f_^^^ (m) for the superposition 
of independent renewal processes with Erlang (k=2) interval distribu- 
tions. The calculations for p = 2,3,4, where p is the number of 
processes superimposed, are exact results using (6.4), (6.7) and (6.9) 
respectively. The result for p = 5 was obtained from the general 
results of Section 5. 

Note that the initial points of the spectra seem to be decreasing 



monotonically with p to the value 



(kTT) 



■, while generally for the 



other 0 ) the spectrum seems to be converging toward the flat spectrum 
of the Poisson limit, i.e. — . 

7T 
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7 . Discussion of Results . 

The results (6.4), (6.7) and (6.8) giving a simple structure to 
the spectra of the superposed processes for p = 2 and k = 2,3 and 
4 are at first quite surprising. However, some such structure could 
be anticipated from (5.15) which gives the spectrum f_^^\o)) as a 
rational function in cos m. It has not been possible, as yet, to 
relate the parameters p and k to the order of the moving average 
and autoregressive structure, or to give simple expressions for the 
spectrum or serial correlations. Knowing the order of the structure 
of the spectrum would be a help in knowing how many serial correlations 
one would have to compute from (5.15) for specific p and k. The 
correlations are computed by calculating f_|_(o)) at a sufficiently fine 
grid and then inverting the series using a FFT algorithm. (See Cooley, 
Lewis, Welch, 1970.) 

There are other theoretical questions raised by the results which 
will be discussed elsewhere. Briefly though, Downton (1971) discussed 
the possibility of looking at this particular superposition process as 
a Markovian progression through pk parallel channels, each channel 
having k stages. This may explain the structure. For p = k = 2 it 
is in fact possible to show that the intervals are not only uncorrelated 
but independent for lags greater than 2. 

We discuss now the computational results. 

In Figure 1 the spectrum is given for the superposition of Erlang 
(k=2) processes for several values of p, the number processes super- 
posed. This was discussed in the previous section. A similar graph 
is given for k = 3 in Figure 2. Again as p increases the initial 
point of the spectrum decrease is converging slowly to the value 1 /tt 
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but the remainder of the spectrum which would obtain for all (x) for 
the limiting Poisson process. 

Figure 3 by contrast has, for fixed number p = 2 of superposed 
processes, graphs for different values k. As k increases and the 
component processes become more regular with C^(X) decreasing, there 
is a larger peak at tt, or equivalently at a period of T = 2. Thus 
the alternation of intervals induced by the superposition becomes 
pronounced in the spectrum of intervals. 

In Figure 4 the case p = 3 is examined for different values 
of k. The tripling due to the superposition produces a more and more 
marked peak at u) = 2tt/3 (period T = 3) as the component processes 
become more regular. 

Formulae such as (5.15) can be developed for non-identical component 
processes. In Figures 5, 6 and 7 several cases of the pooling of a 
Poisson and Erlangian process are given; the marginal distributions 
were given by Lawrance (1973). These are important as models for the 
case in which there is one dominant, fairly regular, process corrupting 
the superposition of many sparse component processes which sum (almost) 
to a Poisson process. Note carefully the vertical scale; departure 
from a flat spectrum is small and would be hard to detect in real data 
unless the Gamma process were very regular. The spectrum of counts would 
probably be much more informative in this case. 

A short table of computed coefficients of variation of the 
marginal intervals X in the superposed process is given in Table 1. 

These came out of the computations of the spectrum. 
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POOLED ERLANG PROCESSES 
Coefficient of Variation 



P =2 3 4 5 6 7 8 9 



.791 


.833 


.860 


.879 


.893 


.904 


.913 


.920 


.714 


.779 


.819 


.847 


.867 


.882 


.894 




.675 


.754 


.802 


.834 










.652 


.740 














.637 


.732 














.627 
















.620 
















.615 

















































TABLE 1 
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8 . Conclusion . 

The fact (Cox and Lewis, 1966, Ch. 8) that the variance-time 
curve and spectrum of counts in the superposition of identical processes 
is a scaled version of the corresponding components is very useful if 
p is known. The crux of most analysis though is that p is not known 
and must be determined. The spectra of intervals computed here should 
be a very useful tool in helping to determine p, and in verifying 
assumptions such as those made in Cox and Lewis (1966, Ch. 8) in 
analyzing nerve-pulse data. 
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CAPTIONS 

Figure 1. Computed spectra of intervals for superposed Erlang renewal 
processes (k=2) for several values of p, where p is the number 
of independent processes superposed. For p = 2, 3, and 4 these are 
exact results using (6.4), (6.7) and (6.9). 

Figure 2. Computed spectra of intervals for superposed Erlang renewal 
processes (k=3) for several values of p, where p is the number 
of independent processes superposed. For all p the results are 
obtained from the computational formulae of Section 5. 

Figure 3. Here the number p of superposed processes is held fixed at 

2, while the parameter k in the pooled Erlang processes is increased 

As k increases the intervals in the individual processes become more 

regular and the peak in the spectrum of the superposition process at 

0 ) = 7T becomes more pronounced. 

Figure 4. Here the number p of superposed processes is held fixed at 

3, while the parameter k in the pooled Erlang processes is increased 

As k increases the intervals in the individual processes become more 

regular and the peak in the spectrum of the superposition process at 

0 ) = 27t/ 3 becomes more pronounced. 

Figure 5. The spectrum of intervals of superposed Poisson and Erlang 
renewal processes is plotted for increasing values of the parameter k 
in the Erlang process. The Poisson and Erlang process mean values are 
held constant. 

Figure 6. The spectrum of intervals of superposed Poisson and Erlang 
renewal processes is plotted for increasing values of the parameter k 
in the Erlang process. The Poisson and Erlang process mean values are 
held constant. 

Figure 7. The spectrum of intervals of superposed Poisson and Erlang 
renewal processes is plotted for increasing values of the parameter k 
in the Erlang process. The Poisson and Erlang process mean values are 
held constant. 
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